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The rapidly growing energy market constantly discovers new alternatives for generating environmen- 
tally-friendly electricity from sustainable energy sources. An externally-heated traveling wave thermo 
acoustic Stirling heat engine is a leading candidate, operating using a variety of viable heat sources. 
Although this engine holds great promise for a low-cost and maintenance-free solution, its current 
reported efficiency still inhibits its market penetration, probably owing to the friction introduced by 
the single moving element - the linear alternator. 

This research quantifies the main parameters affecting the complex thermal-acoustical-electrical 
system, clarifying the critical role of frictional losses. An analytical model has been developed, enabling 
examination of the influence of the critical physical parameters on the electro-acoustic conversion 
efficiency. A measurement method for precise determination of the mechanical friction constant has been 
developed, which enables measuring the friction at the engine’s working frequency. A direct measure¬ 
ment of an engine’s transfer impedance at room temperature enables to find the exact natural frequency 
before field operation, and calibrate external hardware accordingly. A detailed simulation using 
DeltaEC™ indicates that the tight seal gap between the moving piston and its cylinder has a significant 
impact on the system’s overall efficiency. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Efficient Thermo Acoustic Stirling heat Engines (TASE), allowing 
for reliable operation without moving parts or sliding seals, have 
recently stirred a revival of interest. These environmentally- 
friendly, externally-heated engines may operate using a variety 
of viable heat sources, including natural gas, biofuels, open fire or 
solar radiation. Additionally, the lack of need for complex fabrica¬ 
tion processes in the simple TASE structure makes it an excellent 
candidate for commercial solar electricity production. Nowadays, 
the rapidly growing solar-thermal market finds great promise in 
new methods for converting Concentrated Solar Power (CSP) into 
electricity. 

The “traveling wave” thermo-acoustic engine employs the 
inherently reversible Stirling cycle, theoretically capable of 
attaining Carnot efficiency. A practical demonstration of a ther¬ 
moacoustic engine has produced a thermal efficiency (heat to 
acoustic power) of 30% [1] at 725 °C, comparable to that of modern 
internal combustion engines [2] (25-40% heat to mechanical 
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power). Converting acoustic power into mechanical and/or electri¬ 
cal power is no small challenge, and forms the subject of this work. 
A commercial thermoacoustic system has produced over 17 kW of 
acoustic power [3] which was then used to liquefy natural gas 
(Methane liquefies at 115 K at atmospheric pressure) at a rate as 
high as 140 gpd. In 2004, a remarkable thermoacoustic engine, 
converting heat into electricity, was developed by Northrop Grum¬ 
man Space and Technology group with Los Alamos National Labo¬ 
ratory under a NASA contract [4], This engine produced 39 W of 
electrical power at 18% overall conversion efficiency (heat to 
electricity), where a pair of dual opposed linear alternators were 
driven by acoustic power [5], Using Inconel 625 for the hot heat 
exchanger has allowed working temperatures as high as 650 °C 
which enabled the first implementation of a traveling wave 
thermoacoustic engine for efficient electricity production. 

Higher power levels have been demonstrated in 2005 by Luo 
et al. [6], as the nonlinear dissipation losses were reduced in the 
acoustic resonator by employing a tapered resonance tube. In 
2008 [7] more than 100 W electrical output power were produced 
through proper coupling between the linear alternator and the 
thermoacoustic engine, and in 2011 [8] the electrical output power 
was further increased to as much as 481 W with thermal to 
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Nomenclature 


A cross section area, m 2 

B magnetic field, T 

(Bl) alternator constant, N/A 

C compliance, m 3 /Pa 

C e capacitance, F 

F force, N 

/ frequency, Hz 

I, I, total (real), oscillatory part (complex), of electrical 
current, A 

I a acoustic intensity, W/m 2 

I<, Kgas-spring, spring constant (mechanical), gas spring, 
effective, N/m 
L inertance, kg/m 4 

L e electrical inductance, H 

/ characteristic or effective length, m 

m mass, kg 

p, Pm, Pi total, mean part (real), oscillatory part (complex) of 
pressure Pa, (N/m 2 ) 
r Power, W 

R, R e , Rl electrical resistance (internal, external load), Q 
R m mechanical damping, N-s/m 

R r acoustic flow resistor, N-s-m 3 

r space vector, m 

T, T m , T, total, mean part (real), oscillatory part (complex) of 

temperature, K 
t time, s 

U, Ui total (real), oscillatory part (complex), of volumetric 

flow rate, m 3 /s 

u, Ui total (real), oscillatory part (complex), x component of 
velocity, m/s 
V volume, m 3 

Vemf, V h Pint voltage, V 
v velocity vector, m/s 

x axial coordinate, m 

y, z transverse coordinate, m 

Z a , Z alt acoustic impedance (total, acoustic representation of 
alternator impedance), N-s/m 5 

Ze.int, ■Ze.ext, ^e.totai electrical (internal, external, total) impedance, 
Q 


Z e m ,Z e a electrical equivalent of mechanical, acoustic impedance, 

<> 

Z c transfer impedance, A-s/m 

Z m , Z ma mechanical impedance (total, mechanical-acoustical), 
N-s/m 

y c p /c v adiabatic ratio 

ij efficiency 

£ displacement, m 

p, Pm, Pi total, mean part (real), oscillatory part (complex) of 
density kg/m 3 

<t> Pu Ui, <t>A Pu Ui P hase between pi or Ap, and U, 

<d angular frequency 1/s 

Supplementary notation 
9?{} real part 

3{} imaginary part 

|| absolute value, Norm 

()* complex conjugate 

i yrr 

() time average 

A difference 

Subscript 

AE acoustic-electrical 

ALT alternator 

a acoustic 

c characteristic 

CS compression space 

EMF electromotive force 

e electric 

ext external 

EA electro-acoustic 

int internal 

m mechanical 

m, a mechano-acoustic 

REG regenerator 

TA thermoacoustic 

TASE Thermo Acoustic Stirling Engine 

1 oscillatory variable 


electricity efficiency of 12.65% at 650 °C, using 3.5 MPa pressurized 
helium at typical engine working conditions. 

Additional recent projects include the work done at the Energy 
Research Center of the Netherlands (ECN) by Tijani and co-workers 
[9-11) ; the FP7 cooperative project on thermoacoustic technology 
for energy applications led by Spoelstra, Thermeau, Jaworski and 
de Blok [12-14]; the power converter built by Telesz at Georgia 
Institute of Technology [15]; the large-scale-engine research per¬ 
formed by Qiu et al. at the Institute of Refrigeration and Cryogen¬ 
ics, Zhejiang University, Hangzhou, China [16];the investigation of 
low cost regenerator materials conducted by Yu, Jaworski and 
Abduljalil at the School of Mechanical, Aerospace and Civil Engi¬ 
neering, University of Manchester, England [17,18]; and the 
Score-Stove™ project aiming to generate electricity in developing 
countries using thermo-acoustics powered by burning wood. This 
project initiated a broad collaboration between the University of 
Nottingham, University of Manchester, Queen Mary University of 
London (QMUL), Imperial College London and Los Alamos National 
Laboratory (LANL) [19-23] ; Recently, a solar-operated thermoa¬ 
coustic Stirling engine has been reported for the first time to 
produce as much as 200 W of electrical power with overall effi¬ 
ciency of 12.7% (heat to electricity) at 750 °C [24]. The additional 
option of heating the engine with solar-thermal energy has led to 


a compromise of the performance of this prototype compared to 
earlier versions [8], 

Despite its conceptual simplicity and high reliability, nowadays 
further efficiency improvements are necessary for a complete 
TASE-based system to gain a market share. Other well-established 
heat-to-electricity conversion systems have been around for a long 
time, and competing with them, mainly for efficiency, represents 
an important challenge for the TASE system. 

The main objective of this research has been to provide a qual¬ 
itative understanding as well as quantitative analysis of the key 
parameters affecting the complex thermal-acoustical-electrical 
system. A fairly involved trade-off exists between different param¬ 
eters affecting the performance of the system. In particular, we 
address the effect of mechanical friction, known to have a deleteri¬ 
ous impact on the system performance. To better understand this 
effect, an analytical model is developed in this paper. Its results 
will be compared to a detailed thermoacoustic simulation using 
DeltaEC™ in a separate paper. An alternative method of measure¬ 
ment for precise determination of the mechanical friction constant 
was developed and is described. 

Section 2 gives a short background of thermo-acoustic-electrical 
transduction. An analytical model of the involved processes is 
presented in Section 3. Section 4 describes the experimental setup 
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and main experimental results, followed by a short summary and 
main conclusion in Section 5. 

2. Thermo-acoustic-electrical transduction basics 

Thermoacoustic oscillations are the outcome of interactions 
between heat and sound. The development of thermoacoustic 
technology has progressed with many offshoots and branches. 
The first commercial product based on the thermoacoustic effect 
became available in the market approximately 200 years after its 
recognition. In parallel, Stirling cycle engines were investigated 
for almost 150 years, until they found their use in solar powered 
systems, demonstrating the current world record for solar to grid 
efficiency-32% [25], exceeding the previous one - 31.25% [26], 
which can be compared to the highest theoretical efficiency of a 
single p-n junction solar cell - 33.7% [27], 

The principal idea in a TASE is based on a traveling wave passing 
through the gas in a regenerative heat exchanger (regenerator) in a 
sequence of processes: displacement toward a high (heat source) 
temperature, depressurization, displacement toward a low (heat 
sink) temperature, and pressurization. The gas experiences thermal 
expansion during the displacement toward the higher temperature 
and thermal contraction during the displacement toward the lower 
temperature, hence the acoustic power is amplified during the 
cycle; this is the origin of the available-work of the traveling-wave 
engine. The overall efficiency )/ for electricity generation from heat 
depends on the thermoacoustic efficiency i) TA as well as the electro 
acoustic transduction efficiency through their product 

>1 = 1ta'1ea ( 2 - 1 ) 

Wakeland 28] discussed the complexity in reaching the highest 
overall efficiency in thermoacoustic systems by correct impedance 
matching, while Paek et al. [29] as well as Yu et al. [22] suggested 
several experimental verification techniques. In most cases, the 
optimal design is found through some compromise between maxi¬ 
mum output power and overall efficiency. Usually, some degree of 
freedom exists in the integration of a given alternator and a 
thermoacoustic engine. 

The notations used herein below follow closely the ones defined 
by Swift [30,31], A lossless acoustic wave having an angular 
frequency to, related to the primary frequency / through 

co = 2nf (2.2) 

may be described by the following two governing equations - the 
continuity equation 

|° = -V-(pv) (2.3) 

and the inviscid momentum equation 

p(5 + (v-V)v) =-Vp (2.4) 

Using the “Rott lossless approximation” [32-34] means that all 
time-dependent variables are small compared to their mean value 
and only a single frequency is considered (no higher harmonics): 

p(r,t)=p m + 'JJ{p 1 (x)e i “ t } (2.5) 

p(r,t) = p m (x) + 5R{p 1 (r)e te } (2.6) 

r(f,t) = T m (x) + 91{T 1 (r)e te } (2.7) 

U(x,t) = 9?{[/i(x)e tot } (2.8) 

In this approximation the mean pressure p m has no spatial 
distribution whereas the pressure gradients of the oscillating part 
P\ are allowed only in the x direction. The mean temperature T m , 


as well as the mean density p m , are assumed constant in each cross 
section while only the oscillatory parts Ti and p^ can vary in the y 
and z directions. The volumetric flow rate U has only an oscillatory 
part U i since no large mean flow rate is expected. The time 
independent portion of U, referred to as streaming, is only second 
order small, hence would not be discussed here. 

Using this approximation for an ideal gas, the continuity 
equation yields [30] 


P i = - 


yPm 

iaoAAx 


AU, = —AU, 

co C 


(2.9) 


and the momentum equation results in 

Ap, = -irn^^U, = -icoLUi (2.10) 

/I 


while the pressure is related to the volumetric flow rate through 
p 1= -Z a U, (2.11) 

with the definitions of compliance C and inertance L accordingly: 
V 

Compliance - C =- (2.12) 

yp m 


Inertance - L = (2.13) 

Considering a given geometry, the acoustic impedance Z a could be 
specified using combinations of compliances, inertances and resis¬ 
tive elements, as will be explained next. These relations become 
very useful in using the analogy to a lumped electrical circuit. 

An electricity-producing thermoacoustic engine may generally 
be presented according to Fig. 1(a) [4,5], The acoustic elements 
are the compliance C, compression space CS, the inertance L and 
the regenerator REG. Fig. 1(b) describes the elements of the 
alternator, and its connection to the TASE loop. The two pistons 
are each driven by a magnet and inductor, and held by a spring. 
A LVDT is shown for measurement purposes only - to determine 
the location of the piston during operation. 


3. Analytical model 

Fig. 2 describes a simplified acoustical representation for the 
system of Fig. 1 in terms of an electrical analog, where the 
regenerator REG may be represented as a combination of an acous¬ 
tic flow amplifier and a flow resistor R r [35]. The acoustic represen¬ 
tation of the alternator impedance Z ALT includes a contribution 
from two terms: the alternator mechanical impedance as well as 
electrical impedance determined by the alternator’s configuration 
and its main elements. Normally the electrical portion is affected 
by external hardware such as resistive load and tuning capacitor, 
as well as the alternator’s internal structure elements - the mag¬ 
netic field strength, the winding-coil length and diameter, the flux 
closure, and its relative to its equilibrium position. The mechanical 
effect is usually determined by the piston’s cross sectional area, its 
moving mass, the mechanical friction and the alternator spring- 
constant. The exact form and all related details will be elaborated 
later in this paper. 

Three perspectives may be used to analyze the complex system 
response, each having its pros and cons: Fig. 2 shows these per¬ 
spectives, described next, as colored arrows. The TA engine analysis 
mode includes all physical phenomena though it typically requires 
a very complex analysis for a complete inspection; the speaker 
analysis mode is directly available using a simple set of experimen¬ 
tal test procedures though it uses the system in (an almost) 
opposite course of operation; while the EA analysis mode may be 
used to draw exact conclusions regarding the system performance; 
nevertheless, it suffers some ambiguity concerning the incoming 
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Fig. 1 . Schematic representation of a Thermo-Acoustic Stirling engine: (a) complete TASE and (b) alternator in its pressure vessel. 



(b) 


cs 


Pu 


Speaker mode 


Pi.b 


Fig. 2. TASE acoustic representation [35], (a) Complete TA system and (b) alternator 
alone, used as dummy load. 


acoustic pulse characteristics, especially the pressure amplitude, 
the volumetric flow rate and their respective phase differences, 
since the operation generally depends on specific coupling of the 
engine and its load. 

(1) In the TA engine analysis mode (orange), the thermoacoustic 
engine is driven as an engine converting heat into electrical 
power. It is coupled to an acoustic load represented by the 
linear alternator along with its electrical wiring and external 
load. This point of view is quite complex and typically 
numerical tools such as DeltaEC™ or Sage™ are preferred 
over the analytical perspective due to the challenges found 
in capturing the complete picture. 

(2) In the speaker analysis mode (blue) the linear alternator is 
driven by an external AC power supply to produce acoustic 
power. It has an acoustic load defined by the engine geome¬ 
try and acoustic elements. In another form of the speaker 
mode the alternator is separated from the TA engine and 
its connections to the TA engine are blocked. We refer to 
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the alternator operated in this mode as “dummy” load. 
Using this point of view allows for the assessment of each 
acoustic building block. 

(3) In the EA analysis mode (green) there is an assumed acoustic 
pulse characterized by Ap, and U , and their phase difference 
(/) Ap| (J| . In this point of view the electro-acoustic efficiency 
Pea may be estimated. The alternator physical neighborhood 
in its pressure vessel is complex: The piston front face (in the 
compression space) is acted upon by the acoustic wave 
generated by the TASE loop. The moving mass is connected 
with springs to the housing vessel and separated from its 
cylinder with a small seal gap. The back space geometry is 
complex due to the inductor, springs bolts and fixtures, 
and effectively operates as a "gas spring” due to the gas 
volume compressibility. 

Normally, when the system is operated in an engine mode the 
electrical impedance will include an external resistor as well as a 
coil with internal Ohmic resistance and an inductor, connected in 
series to an accompanying power factor-correcting capacitor, 
whereas working in a speaker mode, it should include an AC power 
source instead of the external capacitor and the resistive load. 

Fig. 3 shows simplified schematics of a reciprocating piston 
inside its vessel attached to the TASE. The moving magnet having 
a mass m and piston with cross sectional area A is driven by the 
acoustic pulse, while being mechanically connected to a spring 
having a constant K, and electrically connected to an external load. 
The amplitude and phase of the acoustic pulse Api is simply the 
difference between the pressure oscillations on the front p,/ and 
the back side p, 6 . 

APi =Pi/-Pi,ir (3.1) 

Since the fundamental electro-acoustic transduction analysis is 
complex, it may be useful to begin the discussion with the most 
basic models of the acoustic portion and electrical portion alone, 
while later combining all elements into the full picture, and finally 
extending the model to include additional subtle phenomena. 

Notice that the three different perspectives discussed in this 
section, regarding Fig. 2, are presented in the acoustic viewpoint, 
elaborated next. 

3.J. Acoustic viewpoint 

Considering the force balance on a single oscillating piston, a 
linear acoustic analysis is discussed next assuming all variables 
are not frequency-dependent. In this analysis an acoustic point of 
view is taken by considering a constant cross section - A which 


is acted upon by the oscillatory force - Ap,A, where Ap, is the 
pressure difference across the piston. This pressure difference 
(force over the cross section) is imposed on the piston by the TA 
engine through the volumetric flow rate - LA and the local acoustic 
impedance Z a 

Ap, = -U,Z a . (3.2) 

Once current begins to flow in the alternator coil having a length /, a 
Lorentz force - A p ALT A will be formed due to the magnetic induction 
- B 

, 11 El 

A Palt-—^-- (3-3) 

Furthermore, there is also a pressure drop due to internal 
mechanical forces - A p int due to the piston’s mass - m attached 
to a spring with constant - K, and having mechanical damping - 
R m . The internal pressure drop term is best described by the 
mechanical impedance Z m : 

APint = ^2 = ^2 (J^ + ia)m + Rrn'j ■ (3.4) 

Combining Eqs. (3.2)-(3.4) results in an overall acoustic pressure 
difference balance described in the EA analysis mode as: 

0 = Ap, + A Palt + A p int = -Z a U x + ^ + Z m ^ (3.5) 

Fig. 4 shows a diagram of the associated pressure differences 
described above. 

3.2. Electrical viewpoint 

Because of the acoustic oscillations of the alternator piston 
inside the coil, an electromotive force - V EMF is produced due to 
the change in time of the magnetic flux 

u U ' Bl n fit 

VEMF = Jr ■ (3.6) 

Often, (Bl), the product of the magnetic induction B and an effective 
length / is used as an effective alternator constant. 

Considering the coil as a perfect electrical inductor L e having an 
electrical resistance R e , a voltage V int will develop between the coil 
terminals once a current /, flows through. This can be stated in 
terms of internal electrical impedance - Z eint , while an additional 
power factor correcting capacitor C e may be added in series 

V^int = ~hZ e .i nt = —h + ICOf-e + J • (3.7) 



Fig. 3. Physical picture of the alternator, oscillating inside a pressure vessel, connected to an external electrical load. 
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Fig. 4. Acoustic diagram, showing the pressure difference balance at the alternator- 
TA engine interface. 

The useful voltage Vj which can be derived from the system would 
be consumed by the resistive external load R L 


co r 2n,m 1 

Ia = 2^zJ 0 ^ p i e >' *{ u i e } dt = 2 w {Pi ' u \} (3.12) 

where u\ is the complex conjugate of the oscillatory part of the 
velocity, 9?{} denotes taking the real part of a complex number, 
and Pi is the oscillatory pressure amplitude. In terms of the acoustic 
power - P a , the time-averaged energy flux flowing in the positive x 
direction per unit time (measured in Watts) at a specific location - 
x, through the normal cross section (y-z plane) - is in fact the 
space-averaged acoustic intensity: 

Pa(x)= [ = 

] cross Z 

u section 

= i|p,(x)|.|[; 1 (x)|cos(^ i , 1Ji ) (3.13) 

where <t> p ^u, is the (temporal) phase between the volumetric flow 
rate - Eh and the oscillatory pressure p,. 

Considering Eqs. (3.13) and (3.5) the time averaged acoustic 
power flowing in the TA engine compression space, from the 
acoustic loop into the alternator piston, is given by 



TA engine 



Vi = URl- 

(3.8) 

Ze,ext — Ze,int + Ri. 

(3.9) 


Combining Eqs. 3.6, 3.7 and 3.8 results in an overall electric 
potential difference picture, neglecting any inductive or capacitive 
external elements which may be added later in order to improve 
the power factor. 

Vi = V EMF + Vtot = L ' 1 ^ - /, (Re + icoLe + . (3.10) 

Fig. 5 shows the diagram of the voltage balance described above. 


P.[W] =lsR{Ap 1 u;} = -i 




(3.14) 


while the average electrical power flowing from the alternator into 
the external load is found by considering Eq. (3.10) 


P e [W] =isR{V 1 /*} = 1 


Ml 


Mjl/jIJj-Ke |J,| 


(3.15) 


Combining Eqs. (3.14) and (3.15) results in the following expression 
for the electro-acoustic efficiency, clearly demonstrating how the 
acoustic impedance affects the electro acoustic efficiency (see 

Appendix A.l): 


3.3. Inclusive perspective- Linear alternator as electro-mechanical 
transducer 


91 {Z a } 

1 Iea = 




(3.16) 


The combined overall picture consisting of both acoustic and 
electrical parts is shown in Fig. 6 

The electro-acoustic transduction efficiency^ depends on the 
ratio between the electrical output power - P e and the acoustic 
power - P a produced by the TA engine 


Am — 


Pe 

Pa 


(3.11) 


The (time averaged) acoustic energy carried by a sound wave is nor¬ 
mally stated in terms of the acoustic intensity - I a (energy per unit 
time and unit area) as [36]: 



Fig. 5. Electric-potential impedance diagram for the alternator, and the accompa¬ 
nying electrical connections. 


An alternative form for the electro-acoustic transduction efficiency 
ijea Eq. (3.16), focusing on the effect of the electrical load R L , results 
from the definition of the external electrical impedance Z eext (3.9) 

(Appendix A.2): 


^EA 


n | n i \Ze,ext\ 

Re+RL+^r- 


(3.17) 


Although the wish for low friction R m , minimum electrical resis¬ 
tance R e as well as working at electrical resonance might seem obvi¬ 
ous, indicating the best external load value R L , as well as having 
large (Bl) might be essential for a specific application. Moreover, 
the critical effects of mechanical friction and internal electrical 
resistance on the electro-acoustic efficiency do not allow choosing 
the best external load value a priori [37], Considering operating at 
electrical resonance (Z eex t = R e + Re.ext ) with a given ( Bl ) value of 
50 N/A, Fig. 7 shows the effect of different R e values on rj EA at 
R m = 10 Ns/m, while Fig. 8 shows the effect of different R m [Ns/m] 
values on at R e = 1 £2. 

Knowing the exact values of mechanical friction and internal 
electrical resistance allows choosing the (theoretically) best R L 
value in view of Eq. (3.17). However, indicating the exact R m value 
may require some precaution, as will be discussed next. 

Using the thoroughly explained practice of electro-mechanical 
transduction [38,39], the alternator’s efficiency while operating 
in the speaker analysis mode can be found by the same kind of 
linear analysis of the overall acoustic-mechanical-electrical 
impedance of the experimental apparatus [40]. A compacted form 
for this point of view is shown in Fig. 9: 




















































900 


E. Gonen, G. Grossman/Energy Conversion and Management 88 (2014) 894-906 




z 



Fig. 6. Combined electro-acoustic transduction impedance diagram. 


riEA 



Fig. 7. Internal electrical resistance R e effect on electroacoustic transduction- 
efficiency rjM vs. external load resistance R L , at constant mechanical friction factor - 
R m = 10 Ns/m. 


r/EA 



Fig. 8. Mechanical friction R m effect on electroacoustic transduction efficiency - t] u 
vs. external load resistance R L , at constant electrical internal resistance R ( . = 1 Cl.. 

Notice the reciprocal definitions of electroacoustic transduction 
efficiency between the engine mode Eq. (3.11), and ij^ for the 
speaker mode 



Since the energy- flow changes direction while working in the 
speaker mode relative to the TA engine mode, the pressure balance 
described in (3.5) for the TA engine needs to be modified for the 
speaker mode according to 



Fig. 9. Compact impedance diagram of electroacoustic transduction. 

0 = Api + A p ALT - Ap int 

MB1) /Z m \ (3.19) 

-JT- U \tf +Za )■ 

In the latter case, the energy flow direction is from the linear 
alternator into the acoustic load, following some internal pressure 
drop as well as mechanical dissipation originating from the 
term in (3.4). In the following section the speaker mode will be used 
for evaluating the mechanical friction R m from the experimentally 
measured amplitudes of the current and piston velocity. 

Surely, there are various alternative forms for representing the 
overall impedance diagram in both acoustic framework or 
electrical framework [28,29,22,39,40]; however in this paper the 
following notation is used when considering Figs. 9 and 10: 
Electrical equivalent of mechanical impedance 



Fig. 10. Overall electrical representation . 
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- m 2 

m 7 


Electrical equivalent of acoustic impedance 



Internal electrical impedance 

Ze.int — Re + i(OL e + 

lC0C e 


Overall impedance 

7 7 1 

^e,total — ^e,int "r -j j 

75- + 7r 


— J?e H - icoL e ———h 


(Blf 


R m + itom + £ + A Z 0 


(3.20) 


(3.21) 


(3.9a) 


(3.22) 



Fig. 11. Qdrive's 1S102M/A linear reciprocating alternator. 


Yet another simplification may be introduced by considering the 
mechanical-acoustical impedance Z m a 

Z m ,a = Z m + A 2 Z a . (3.23) 

Z e ,to t a,=Z eiiBt +^-. (3.24) 

Considering the above mentioned shortcuts, the current in the loop 
h is related to the externally measured (or supplied) voltage V^, 

Fig. 10. 

V, =hR L . (3.25) 


different environments. The STAR has a complex design, having 
non-axi-symmetric “interdigitating interfaces” between moving 
and fixed magnetic circuit elements [42], Nevertheless, its key 
operational principles are common to all linear alternators, where 
a conversion of power occurs from an oscillatory mechanical 
energy into an electrical energy in the form of alternating current. 
Similar to common linear alternators [43], this device has a Vi kg 
piston reciprocating back and forth on a linear trajectory, inducing 
a change of magnetic flux in the surrounding inductor wire loop. 
The voltage across the loop is proportional to the rate of change 
of the magnetic flux through the wire loop according to Faraday’s 
law. 


URl + UZeioiai — 0. (3.26) 

The complexity in achieving the highest overall efficiency in a ther¬ 
moacoustic system by correct impedance matching has been dis¬ 
cussed thoroughly by Wakeland 28], while several experimental 
verification techniques have been suggested by Liu and Garrett 
[37], Paek et al. [29] and Yu et al. [21,22] . Under most cases the 
optimal design is found through some compromise between maxi¬ 
mum output power and overall efficiency. Moreover, usually there 
is some degree of freedom in the integration of a given alternator 
and a thermoacoustic engine. Recent publications disclose two 
noteworthy coupling methods of electro-acoustic coupling 7,10], 
Furthermore, it can be found that the volumetric flow rate U, is 
related to the input voltage Vt through (Appendix A.3): 


Z eM (Z m +A 2 Z a ) + (Bl ) 2 

Based on these analytical considerations, proper selection of imped¬ 
ance values is crucial. Some degree of freedom remains available for 
optimization of a complete electricity-producing thermoacoustic- 
system, despite technological challenges such as minimizing fric¬ 
tion, maintaining tight seal gaps and leaks prevention, as well as 
practical limitations related to piston stroke, mechanical springs 
elongation and dead volumes. However, accurate evaluation of crit¬ 
ical impedance values is essential for comprehensive optimization. 

4. Experimental setup and results 

This section starts with a short description of a simple prelimin¬ 
ary test procedure for evaluation of mechanical friction, followed 
by a general description of our experimental apparatus, concluding 
with detailed impedance assessment and validation of the analyt¬ 
ical model. 

During the experimental stage a custom design of STAR 
1S102M/A linear alternator [41], Fig. 11, was used in several 


4.1. Preliminary testing 

The simplest approach for mechanical friction determination is 
by following the free decay of a linear alternator’s movement [44], 
placed in standard atmospheric conditions and held firm by its 
stand. It would seem rational to implement this by giving the alter¬ 
nator a jolt, either mechanical or electrical, thus disturbing its 
equilibrium and causing the piston to oscillate at its natural fre¬ 
quency until brought back to rest by the friction. However, this is 
impractical when the alternator is mounted in its pressure vessel 
and the TASE filled with gas. Therefore, to apply the Resonance 
Decay method, the alternator’s natural resonance frequency must 
be identified first, for example by locating the peak of a frequency 
response measurement, as described next in Section 4.2; while 
later the piston is driven at its resonance frequency using an exter¬ 
nal power supply, until the measurement process begins as the 
driving power source is shut abruptly. Fig. 12 shows the theoretical 
approximate displacement fit (red ) from Eq. (4.1) [44], for the 
experimental damped harmonic oscillations (blue), indicating 
mechanical friction characterized by R m = 6.9 [NsitT 1 ], measured 
in the laboratory at 1 atmosphere of air at room temperature. 

£jit(t) = |£| sin(cot)e- Rmt/2m . (4.1) 

Note that this preliminary test method does not provide the full fre¬ 
quency dependence of friction. In a practical system it is necessary 
to determine the friction at the operating frequency, which may be 
different from the natural frequency. Yet, in our preliminary phase 
experiments, performed approximately one year after the alternator 
was purchased, the measured value was somewhat different than 
the manufacturer’s specified value of R m = 5.5 [NsitT 1 ], which 
seems reasonable since some inevitable wear has meanwhile 


1 For interpretation of color in Figs. 7,8,12, and 14-17, the reader is referred to the 
web version of this article. 
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Fig. 12. Alternator piston’s free decay at 39.8 Hz. 


occurred. 


4.2. Impedance assessment and validation testing 

To confirm the validity of the complete analytical model 
described above, as well as assessing critical impedance values, a 
set of experiments was conducted. The experimental apparatus 
details are presented briefly, followed by a specific description of 
the laboratory setup, closing with highlights of the major results. 

Our experimental system shown in Fig. 13 consisted of a TASE 
equipped with an alternator, as described in Fig. 1. 

During the experimental testing phase, the system was driven 
both as an engine, Fig. 1(a), converting heat into electrical power, 
or as a speaker, producing acoustic power driven by an electrical 
AC power supply. For that purpose, an adaptable experimental 
system was fabricated, assembled, and operated at different condi¬ 
tions. The versatile system allowed for a good comparison between 
sequential sets of experimental testing. In an additional “dummy” 
speaker configuration, the compression space was blocked, thereby 
separating the alternator section from the TASE loop, hence modi¬ 
fying the compression space acoustic impedance, as well as 
detaching the regenerator and heat exchangers from the available 
acoustic power. 

In the speaker mode, the system typically consisted of a 
dual-opposed linear alternators driven by an external programma¬ 
ble AC Source, Chroma model 61504 [45], attached to either a 
thermoacoustic Stirling engine, or to a simplified “dummy" load 


volume. The entire experimental apparatus included two types of 
electricity generating ThermoAcoustic Stirling Engine models, 
termed TASE-200 and TASE-B. Each TASE may be operated as an 
engine - converting heat into electricity, or in a reverse mode of 
operation - employing electricity for heat pumping under a tem¬ 
perature difference at the desired working conditions. During the 
engine mode of operation, the three heat exchangers are used to 
supply and remove heat from the engine, as the thermoacoustic 
oscillations are driving the alternators producing electricity, 
consumed by an external electrical load. In the reverse mode of 
operation a pair of dual opposed linear alternators are powered 
by the external electrical power supply, producing acoustic power 
which may be used to pump heat from the low temperature to the 
high temperature heat exchanger. These models of thermoacoustic 
engines operate much like earlier versions described in the litera¬ 
ture [4,5,15], During the electro acoustic transduction testing, the 
alternators were connected in series to the external power supply, 
delivering the AC current at predetermined frequencies, as the 
external voltage across the alternators was kept constant through¬ 
out the frequency sweep. 

Two types of pressure sensors were installed along the loop at 
critical points, comprising Endevco 8510B [46] piezoresistive pres¬ 
sure transducer, as well as a PCB Piezotronics [47] piezoelectric 
subminiature ICP® pressure sensor model 1050 2 , along with 20 
thermocouples situated in the thermodynamic section, spanning 
the heat exchangers, regenerator and the thermal buffer tube 
walls. The piston’s displacement was measured using LVDT Macro- 
sensors model CD-375-250 [48] while the current and voltage 
readings were fed directly into the NI compact RIO system [49], 
All the required data acquisition and preliminary analysis were 
performed in real time using the Labview™ program. During each 
test, the alternator piston was excited at various driving frequen¬ 
cies as the supplied voltage remained constant. The calculated 
velocity amplitude found by measuring the displacement ampli¬ 
tude at each driving frequency, may be compared with the velocity 
u deduced from the analytical model (Eq. (4.2) derived from 
Eq. (3.27)) for a known piston area, as shown in Fig. 14. 


(BQV i 

Z eM (Z m +A 2 Z a ) + (Bl) 2 ' 


(4.2) 


To calculate the velocity iq the following procedure was used: the 
density p m was derived from the helium properties at fill pressure 
p m of 30 bar at room temperature: the parameters K, m and (Bl) 




Hot heat^xchanger 


Thermal 
buffer tube 
Linear 
alternators 


Regenerator 
heat 
exchangers 

Compliance 


Inertance 


(a) 


(b) 


Fig. 13. Experimental setup of the tested system, (a) Actual TASE Photograph and (b) schematic system illustration . 
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Fig. 14. Alternator's piston frequency response with 3 MPa helium gas. 


were taken from the alternator’s data sheet; the electrical values IA, 
R e , L e and C e were directly measured using standard tools, and used 
to calculate Z ein t from (3.8) . The acoustic portion Z a required for Eq. 
(4.2) was evaluated considering Fig. 2, Eq. (2.12) and (2.13) using 
the geometrical dimensions of the apparatus, having a 2" diameter 
piston and compression space volume of 5E-5m 3 . The inertance 
had a cross section area of lE-4m 2 , as its lengths for the dummy 
and TASE200 were L dummy = 0.14 m, L TASE2 oo = 0.51 m respectively, 
while the TASE200 included additional effective compliance volume 
Vtase2oo = 5.11E-5m 3 . To evaluate the mechanical friction factor R m 
from this velocity-frequency curve we guessed various values of R m 
and used them to evaluate Z m from (3.4) until a good correlation 
between the calculated and the measured values of the velocity 
amplitude was obtained. The results presented in Fig. 14 qualita¬ 
tively confirm the analytical model, though a slight discrepancy 
may be observed at high frequencies. The left peak is mainly con¬ 
tributed by the electrical portion, while the right peak is mainly 
dominated by the mechanical-acoustic portion. 

Comparing the two techniques described above, the frequency 
response method and the alternator’s free decay method, each 
has its own advantages. The free decay method is more focused 
on mechanical effects, while the frequency response method 
depends also on the acoustic impedance Z a as well as on electrical 
parameters. Moreover, in the free decay method, the effect of fre¬ 
quency-dependent phenomena on the alternator’s impedance is 
concealed. This may become misleading, for example in a thermoa¬ 
coustic engine, where the working frequency is normally different 
from the alternator’s natural frequency. The weakness of the sug¬ 
gested method for evaluating R m by correlating the measured data 
with the velocity calculated from Eq. (4.1) originates from the 
requirement of using some alternative techniques in order to attain 
prior knowledge of the exact values of the electrical parameters - 
R e , L e , C e . 

Already during early testing stages, the frequency response 
method yielded a very high mechanical friction factor R m , about 
three times larger than the free decay method. This discrepancy 
may originate from different oscillating flow patterns in the pis¬ 
ton’s vicinity due to dissimilar acoustic environments between 
the two experiments; therefore a more accurate method was 
sought. Notice the main difference in the experimental procedure, 
where in the open air free decay method the alternator was fixed in 
an open test bench, while during the frequency response set of 
experiments the alternator was surrounded by a helium atmo¬ 
sphere inside a pressure vessel at varying pressures. Nevertheless, 
the two configurations tested had the same alternators, but located 
inside different acoustic environments. 


To resolve the inferred difference in friction factor, as well as 
further validate the analytical model effectiveness, the data was 
evaluated in light of another measurable quantity - the transfer 
impedance Z c , defined here as (Eq. (3.20)): 


7 _ |/| _ Z m ,a 

c - M (Bl) 


(4.3) 


Since the transfer impedance is not affected by the electrical imped¬ 
ance, while enabling simple exploration of different working 
conditions, it may become more adequate for highlighting the 
acoustic-mechanical effects. In order to best estimate the friction 
factor at the working point of a thermoacoustic engine, it was found 
that the transfer impedance method was the preferred method 
among the three options discussed in this paper. 

Several different acoustic environments were experimentally 
attainable, depending mostly on the gas free volume V, as well as 
the mean fill pressure p m , forming together an almost adiabatic 
“gas spring” having an effective spring constant - K ga s-spring acting 
on the piston face at a cross section area A p : 

K ~ yPnAp (4 4'j 

‘'gas-spnng ~ ,, • 


By changing either the mean pressure while keeping the volume of 
the Dummy fixed at V dummy = 3E-5 m 3 (Fig. 15), or by changing the 
free volume while keeping the 30 bar mean pressure (Fig. 16) - 
different peak locations are evident in the transfer impedance Z c 
frequency response curve. These changes are expected considering 
an effective spring constant K e ff affecting the magnitude of Z ma on 
the right hand side of Eq. (4.3): 

K e ff ~ K + Kg as _ spr i n g. (4.5) 

Hence, the spring constant K used in Eq. (3.4) for Z m should be 
replaced by K e fj according to Eq. (4.5). 

The transfer impedance method described above allows for a 
straightforward comparison between the analytical model and 
the laboratory measured values, keeping the mechanical parame¬ 
ters constant, while allowing for several system configurations, 
over a broad frequency range as well as different fill pressures. 
Moreover, by measuring an engine’s transfer impedance at room 
temperature one can predict the operating frequency before its 
usage in the field. This may become useful for fine tuning of exter¬ 
nal hardware such as balancer for vibration cancellation, capacitor 
and starter. Yet another possible benefit from the transfer imped¬ 
ance measurement is an estimation of an effective (Bl) product as 
well as a measurement of L e at typical working conditions (1, co, 
{), as will be demonstrated in a following paper. 

The correlations between the analytically predicted values 
(solid lines) and the laboratory measured results (experimental 



Fig. 15. Mean pressure effect on the measured transfer impedance Z c , demonstrat¬ 
ing different frequency curves due to altered stiffness (gas spring effect). 
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1/Zc [m A 1 s’ 1 ] 



Fig. 16. Mutual fit (solid lines) for the measured transfer impedances Z c of TASE 200 
(+) and dummy load (*). 


1/Z C [A - 1 s- 1 m] 



Fig. 17. Simulated effect of different gap heights (10-90 pm) on transfer imped¬ 
ance Z c . 


values) in Fig. 16 seem reasonable for both system configurations - 
TASE200 (+, red) and dummy (*, green). However, the mechanical 
friction inferred from the fit shown in Fig. 16, R m = 45 [Ns/m], is 
higher by a factor of ~6 than the one determined in the natural free 
decay. Furthermore, the respective curve fitting for both configura¬ 
tions was unsatisfactory despite numerous parametric curve fitting 
trials, suggesting that a frequency dependent effect might be 
missing in the analytical model. It might be possible that viscous 
resistance as well as blow - by in the tight seal gap are the 
dominant frequency dependent phenomena to be considered. 

This has initiated the process of extending the analytical model, 
to include the more subtle issues which may have been unnoticed 
before, using standard methods. 

To account for the discrepancies, a numerical simulation using 
DeltaEC™ [50] was performed including viscous losses as well as 
blow - by in the seal gap. In the simulation even negligible effects 
such as thermal-relaxation dissipation losses [30] were considered. 
In the simulation, the externally supplied voltage V, is kept con¬ 
stant while the current, the pressure amplitude, and the piston 
velocity are allowed to vary as the frequency is changed. For each 
gap height a frequency sweep was carried out while forcing bound¬ 
ary conditions suitable for a 30 bar helium closed chamber. Fig. 17 
shows the frequency response of the transfer impedance Z c for five 
different gap heights. As the gap gets smaller the seal becomes 
tighter and the gas spring effect becomes more dominant, shifting 


the resonance peak from its mechanical resonance at about 40- 
100 Hz. When the gap is open no pressure difference can remain 
between both piston sides, making the gas spring effect negligible, 
whereas in a tight configuration, the gas spring effect as well as vis¬ 
cous losses inside the gap become dominant. 

Following a fine-tuning of the simulation to correlate with the 
experimental results the analytical model was extended to include 
the additional losses related with the seal gap. These results will be 
described in a following paper. 


5. Summary and conclusions 

One leading solution for a sustainable energy production is the 
traveling wave thermo acoustic Stirling heat engine; yet, its prac¬ 
tical efficiency limitations seem to delay its market penetration. 
Despite its great promise as a highly reliable device, currently no 
such engine is available commercially, as the technology employed 
is not yet mature enough. It appears that one of the main obstacles 
for wide-ranging acceptance is the friction occurring in the only 
moving element - the linear alternator’s piston flexure seals. 

This research has provided a quantitative analytical examina¬ 
tion of the influence of the critical physical parameters on the elec¬ 
tro acoustic transduction efficiency t] EA , revealing the severity of 
frictional effects present in the system. A novel measurement 
method for precise determination of mechanical friction constant 
- R m has been developed, indicating a mismatch in friction factor 
evaluation at the engine’s working frequency, between the two 
methods described above - the free decay and the frequency 
response methods. Moreover, using a frequency independent R m 
does not allow obtaining a satisfactory respective curve fitting 
for the experimental systems studied due to the involved complex 
interactions. Additionally, a direct measurement of an engine’s 
transfer impedance at room temperature enables to find the exact 
natural frequency before field operation, and calibrate external 
hardware accordingly, leading to practical benefits . 

In this study, two experimental systems (TASE200 and 
dummy) have been extensively explored at a variety of fill pres¬ 
sures to gain greater confidence in the correlations between the 
measurements and the analytical model. The two systems were 
tested repeatedly at a variety of frequencies and detailed 
measurements were conducted to evaluate the operating parame¬ 
ters. These experimental results, devised from a robust set of 
experimental testing, were used to calculate the relevant imped¬ 
ances. In addition, a detailed thermoacoustic simulation using 
DeltaEC™ was performed in order to verify the effect of the seal 
gap on the transfer impedance. Good agreement between the 
experiments and simulations was observed, while in a following 
paper we will extend the analytical model to include the effects 
of viscous friction in the tight seal gap, resulting in better match 
to the experimental values. 

Public acceptance of products based on electricity generating 
thermoacoustic engines is expected to grow in the near term, as 
higher efficiencies are expected to be attained through proper 
implementation of specific elements, accessible through precise 
determination of the physical parameters present in the actual 
structure. Improved designs, based on the considerations explored 
in this paper, may lead to enhanced system performance, mostly in 
terms of the product overall cost per Watt, as well as the envisaged 
higher power levels. 
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Appendix A 


A.3. Velocity calculation- speaker mode 


A.l. Electroacoustic efficiency - acoustic impedance effect 


Considering Eqs. (3.23)-(3.25) and (3.26) 


Substituting the acoustical and electrical powers from (3.14) 
and (3.15), respectively, into the definition of the electroacoustic 
transduction efficiency (3.11) yields: 


V, =I l R L = -hZ efotal = -J, (Z eM + g) 

=► U =- 


(A3.1) 


Rea : 


f X{U,I]}-R e \h\ 2 


From Eq. (3.5): 


(Blfll 

a u : 



i u ; 

}- 

. p Ll II 
(A U\ 

(Bl) 

A 

3 J { 

ni 

<4J 

\+¥ 


7 Zm 
° ~A* 


(Al.l) 


(A1.2) 


and 


Recalling Eq. (3.10) and changing variables from the volumetric 
flow rate Uj to velocity iq yields: 

= V EMF + V tat = u, (Bl) - /,Z, int (A3.2) 

Substituting Eq. (A3.1) into (A3.2): 


(Bl)ut = l x Z ti tat + V, = hZ eMt - h(Z e , iBt + g) = - , ^ r 

=>iq = -|® 

Z-m.ci 


(A3.3) 


Re^-hr^Re 
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A 2 

W 
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A 2 / (Bl) 

7 -m' 2 


Z„-- 


Substituting in Eq. (Al.l): 


Rea : 


H z °-*y- R 'w i z °- z #i 
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*{Z„} + f 

and considering the expression for Z m , Eq. (3.4): 

i2 


Rea : 


®{z«)-)- R w 


«{Z 0 } 


(A1.3) 


(A1.4) 


(A1.5) 


A.2. Electroacoustic efficiency - engine mode 

Considering Eqs. (3.10) and (3.25) while recalling the definition 
for the external electrical impedance Z eext (3.9) 


1/ _ U t (BI) I 7 

Vi = UR l 


Ui (Bl) 


— h(RL H - Z e jnt) — I\Z e ,ext 


(Bl) I, (Bl\ 2 1 

A U, {A) Z e , ext 


(A2.1) 


and substituting into the definition of the electroacoustic transduc¬ 
tion efficiency according to Eq. (Al.l) yields: 



9*] 
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Evaluating the complex terms 
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= 31 
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(A2.2) 


(A2.3) 


and substituting back into Eq. (A2.1) yields: 


R acoustic-electrical 


fte+ftf. _ Re _ 

|^e,exf| l^e.ext I 

Re+RL _i_ Rm 
l^e.extP 



(A2.4) 


and concluding with the voltage V A from (A3.1): 

u V. (Bl) (Bl)V i 

' Z ^ + ( €* Zma At Zm,a + (Bl ) 2 
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